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ABSTRACT 


Mulej,  David  Joseph,  M . S . M . E .,  Purdue  University,  May,  1975* 
The  Velocity  of  the  Interface.  Major  Professor:  V.  W.  Gold- 
schmidt. 

A unique  measure  of  the  lateral  component  of  the  velo- 
city of  the  turbulent-nonturbulent  interface  was  defined  and 
expressed  in  terms  of  the  crossing  frequency  and  intermit- 
tency,  both  parameters  of  the  turbulent  flow  field.  The 
velocity  of  the  interface  was  compared  for  different  turbu- 
lent shear  flows  and  was  found  to  show  a trend  of  similarity. 
Experimental  measurements  of  the  velocity  of  the  interface 
were  taken  in  a two  dimensional  plane  jet  and  found  to  cor- 
respond with  the  predicted  values. 

Other  characteristics  of  the  interface  were  investigated 
and  reported  in  this  work.  It  was  noted  that  the  velocity 
of  the  front  and  back  of  the  interface  were  not  the  same. 

A measure  of  the  amount  of  folding  in  the  interface  as  a 
function  of  location  within  the  jet  was  also  reported.  It 
was  found  experimentally  that  the  amount  of  folding  increases 
directly  with  the  intermittency . 


INTRODUCTION 

This  work  is  a continuation  and  an  outgrowth  of  work 
being  done  in  conditional  measurements  of  turbulent  shear 
flows.  The  work  by  this  author  is  a natural  extension  of 
the  work  by  Pete  Jenkins  (197*0  in  a two  dimensional  plane 
jet.  The  main  investigation  in  this  work  concerns  itself 
with  the  measurement  of  characteristics  of  the  turbulent- 
nonturbulent  interface. 

This  work  will  present  a unique  measure  of  the  velocity 
of  the  interface.  This  will  be  a measure  of  the  velocity  of 
the  interface  in  the  lateral  direction,  which  constitutes  a 
different  view  of  the  interface  than  is  currently  in  the 
literature . 

There  will  also  be  a discussion  of  other  measures  of 
the  interface  including  folding.  The  folding  of  the  inter- 
face, as  discussed  by  Paizis  and  Schwarz  (197*0 » is  defined 
as  the  action  by  which  a region  of  turbulent  fluid  has  non- 
turbulent  fluid  between  it  and  the  fully  turbulent  region. 
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I ANALYSIS 

1 . Review  of  Previous  Work 
The  turbulent  interface  has  been  investigated  for  a 
number  of  years.  Corrsin  (19^3)  is  credited  as  being  the 
first  investigator  to  observe  via  a hot  wire  probe  the 
intermittent  character  of  a turbulent  shear  flow.  Townsend 
(19^8,  19^9)  introduced  the  concept  of  intermittency,  which 
is  defined  as  the  fraction  of  time  that  the  flow  is  turbu- 
lent. He  took  intermittency  measurements  in  the  wake  of  a 
circular  cylinder  and  found  it  to  obey  a similarity  relation- 
ship. 

The  first  major  work  discussing  the  interface  was  that 
of  Corrsin  and  Kistler  (195^)  which  among  other  things 
investigated  the  thickness  of  the  interface.  The  superlayer, 
as  the  interface  was  termed  by  Corrsin  and  Kistler,  was 
found  to  have  a thickness  based  on  similarity  grounds  propor- 

o i 

tional  to  the  Kolmogorov  microscale,  (vv^)4, 

Townsend  (1956)  compiled  much  of  the  early  intermittency 
data  in  free  shear  flows  and  later  (1966)  proposed  a model 
for  the  interface.  His  intention  was  to  propose  a theoreti- 
cal model  for  the  interface  that  would  explain  its  behavior, 
especially  the  relationship  between  the  interface's  highly 
irregular  nature  and  the  entrainment  of  ambient  fluid. 
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Townsend  thereby  proposed  a visco-elastic  model  for  the  mo- 
tion of  the  interface. 

The  technique  of  conditional  averaging  is  a relatively 
new  experimental  method  for  investigating  the  intermittent 
regions  of  shear  flow.  Early  workers  in  this  field  include 
Kibens  ( 1968) , Kovasznay , Kibens , and  Blackwelder  (1970),  and 
Jenkins  (197*0  among  others.  The  basis  for  these  measure- 
ments is  the  ability  to  detect  the  presence  of  the  interface. 
Jenkins  (197*0  measured  the  velocity  and  temperature  profiles 
with  respect  to  the  interface  in  a heated  two  dimensional 
plane  jet.  He  also  presented  an  up-to-date  literature  re- 
view. 


Of  interest  is  the  work  of  Phillips  (1972)  and  that  of 
Paizis  and  Schwarz  (197*0.  Phillips  was  primarily  concerned 
with  the  part  played  by  the  interface  in  the  action  of  the 
entrainment  of  ambient  fluid.  He  proposed  a pseudo -Lagrang- 
ian  description  of  the  entrainment  process.  Paizis  looked  at 
the  shape  of  the  interface  and  noted  the  significant  amount 
of  folding  which  occurs  in  the  interface.  He  also  measured 
the  characteristic  scales  and  convection  velocity  of  the 
interface  for  a turbulent  wall  jet. 

Most  investigators  have  looked  at  the  interface  as  a 
thin  surface  separating  the  turbulent  and  nonturbulent  flows. 
That  surface  is  convected  downstream  by  the  mean  flow.  In 
the  work  that  is  now  presented  the  interface  in  a plane  jet 
is  observed  from  a fixed  position  in  the  shear  flow.  The 
relative  lateral  motion  of  the  interface  is  then  noted.  The 
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motion  is  periodically  inward  and  outward  from  the  point  of 
reference  (without  accounting  for  the  axial  displacement). 
This  view  is  the  basis  for  the  measurement  of  the  (lateral) 
velocity  of  the  interface. 


2 . Derivation  of  Lateral  Velocit 


probe  1 
probe  2 


When  the  probes  are  placed  in  the  intermittent  region 
of  the  jet,  their  output  will  be  similar  to  that  shown  in 
Figure  2.  From  these  types  of  signals  there  are  two  terms 
which  may  be  defined:  the  intermittency,  Y,  and  the  cross- 
ing frequency,  fY. 
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r 

The  intermittency  of  the  flow  is  defined  as 


or  as  the  fraction  of  time  the  flow  is  turbulent.  Figure  3 
shows  an  idealized  output  of  hot  wire  probes  in  the  inter- 
mittent region.  From  Figure  3 the  terms  of  Equation  (l)  are 

clearly  shownt  T is  a large  period  of  time,  t . . is  the  j 

"t  h 

portion  of  time  the  i trace  appears  to  be  turbulent,  and  n 
is  the  number  of  turbulent  bursts  during  the  time  T. 

Using  the  traces  shown  in  Figure  2 it  is  relatively 
obvious  when  the  flow  is  turbulent  and  nonturbulent , and 
also  quite  straight  forward  to  measure  the  intermittency. 
However,  when  the  intermittency  is  measured  electronically, 
the  problem  is  more  difficult.  Any  electronic  method  which 
may  be  used  to  measure  intermittency,  such  as  the  method 
described  by  Jenkins  (197*0*  is  usually  verified  and  cali- 
brated by  using  visicorder  traces  of  the  hot  wire  probes. 

At  present  the  intermittency  of  the  flow  in  a two  di- 
mensional plane  jet  is  well  understood.  Figure  4 shows  the 
data  by  Jenkins  (197*0  in  the  present  setup  and  compares  it 
with  the  data  taken  by  Bradbury  (1965)  and  Heskestad  (1965)* 
The  other  term  which  may  be  deduced  from  Figure  2 or  3 
is  the  crossing  frequency,  sometimes  referred  to  as  the 

I bursting  frequency  when  speaking  of  boundary  layer  flows. 

The  crossing  frequency  is  defined  as 

I fY  = Lim  2 (2) 


Jenkins 
x/D  - -5 

O 

35 

C 

^5 

A 

55 

□ 

Bradbuij 

• 

Heskestad 

■ 
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The  crossing  frequency  has  also  been  measured  in  the 
present  setup  and  this  data  by  Jenkins  (1974)  is  shown  in 
Figure  5- 

The  velocity  of  the  interface  in  the  lateral  direction 
may  be  defined  using  Figure  3 as 

= xf  (3) 

and  at  a point  in  the  flow  field  as 


V. (x,y)  = Lim  ^ 
1 Ay-0  ^ 


which  is  equivalent  to 


dy; 


V.  = — - 
i dt 


where  y^  is  the  lateral  position  of  the  interface. 

If  ^y  is  sufficiently  small  so  that  n^  - n£,  then 


- 

At  = ? - - 1 

2n 


Subatituing  into  (4) 


V.  - Lim 


2nAV 


- MJ  XUJ 

*y~°  4ji.i  - 


j=l 


k=l 


and  using  the  definitions  of  V and  fy 


V\  = Lim 


.y— -0  Y1  ' v2 


V.  = Lim 


2f, 


i ■ 


V.  = 
i 


_2f 

zy 

dV/dy 


(^) 


(5) 


(6) 
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Using  the  data  of  Jenkins  (197*0  (and  repeated  in  part 
by  the  author)  for  Y and  f^ . che  velocity  of  the  interface 
is  found  according  to  Equation  (6)  and  is  shown  in  Figure  6. 
The  velocity  of  the  interface  is  made  dimensionless  by  the 
terms  which  arise  naturally  in  the  analysis.  The  data  of 
Figure  6 is  replotted  in  Figure  7 in  terms  of  Y rather  than 
y/b.  This  gives  a better  view  of  the  distribution  of 
within  the  intermittent  region. 

Various  methods  were  attempted  in  deriving  from  the 
available  data  using  Equation  (6).  The  problem  was  seen  to 
be  accurately  obtaining  the  derivative  of  the  intermittency . 
Initially  a high  order  polynomial  was  used  to  approximate 
the  intermittency  function.  This  agreed  well  in  magnitude, 
but  not  in  the  derivative.  Then  an  exponential  form 

Y = exp( -A (y/b  - c)n) 

was  used  to  fit  the  data  with  A,  c,  and  n determined  through 
curve  matching.  This  choice  gave  easily  obtained  derivatives 
which  appeared  as  a well  behaved  function.  Finally  a graph- 
ical method  was  used  to  extract  the  derivative  from  the 
intermittency  data  directly.  These  results  show  that  the 
exponential  fit  does  not  work  well  for  the  derivative.  Fig- 
ure 8 shows  the  derivative  of  the  intermittency  as  obtained 
from  the  data  of  Figure  4. 


J 
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3 . Comparison  With  Data  From 
Other  Shear  Flows 

Earlier  researchers  have  spoken  of  an  entrainment  vel- 
ocity (the  velocity  at  which  nonturbulent  fluid  is  entrained 
into  the  flow,  Phillips  (1972)),  of  an  interface  convection 
velocity  (the  velocity  of  the  interface  in  the  direction  of 
the  mean  flow,  Bradshaw  (1967)),  and  the  velocity  of  the 
fluid  at  the  interface  (see  Jenkins  (197^)).  Of  concern  now 
is  the  lateral  component  of  the  velocity  of  the  interface, 

V.  . 

l 

The  data  necessary  to  determine  can  be  found  from 
other  shear  flows  measured  by  earlier  investigators.  Inter- 
mittency  and  crossing  frequency  data  have  been  published  by 
Thomas  (1973)  for  the  plane  wake  of  a cylinder  and  by  Kibens 
(1968)  for  a turbulent  boundary  layer.  The  velocity  of  the 
interface  has  been  calculated  from  this  information  using 
Equation  (6)  and  is  compared  in  Figure  9 to  that  found  in 
this  present  work  for  a plane  jet. 

In  order  to  solve  for  the  intermittency  data  for  the 
wake  and  the  boundary  layer  were  first  approximated  by  a 
smooth  function  chosen  to  fit  the  data  (but  without  checking 
for  complete  agreement  with  the  derivative).  The  points 
shown  in  Figure  9 represent  values  taken  from  the  fy  data. 

To  make  a meaningful  comparison,  the  velocity  of  the 
interface  was  made  nond imensional  in  each  of  the  flows  by 
f.m>  the  maximum  crossing  frequency  and  Ay*,  a measure  of 
the  width  of  the  intermittent  region,  given  by 
Ay*  = y( Y=0 . 01 ) - y(Y=0.99) 


(7) 


Figure  9 shows  a plot  of  the  nondimensional  velocity  of 
the  interface  as  a function  of  intermittency.  It  shows  a 
definite  trend  of  similarity  in  the  velocity  of  the  inter- 
face for  different  shear  flows. 

It  should  be  noted  that  although  for  the  boundary 
layer  and  the  plane  wake  were  obtained  as  stated  above,  the 
corresponding  for  the  current  work  was  obtained  by  graph- 
ically differentiating  the  intermittency  data  and  reading  off 
the  values  of  fy.  It  is  felt  that  the  results  shown  in  Fig- 
ure 9 are  not  affected  by  this  difference. 


2.5 


O Plane  Jet  (present  work) 
□ Plane  Wake  (Thomas) 

A Boundary  Layer  (Kibens) 


Figure  9 Velocity  of  Interface  vs.  y for  Different 
Shear  Flows 
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II  EXPERIMENTAL  PROCEDURE 


1.  Experimental  Setup  and  Characteristics 

The  experimental  setup  used  for  this  work  was  essential- 
ly the  same  as  used  by  Jenkins  (1974).  The  two  dimensional 
plane  jet  (see  Figure  10)  was  powered  by  a 0 . 56  kW  (0.75  bp) 
blower.  The  jet  has  the  capability  of  heating  the  flow  by 
use  of  a 4.0  kW  mesh  wire  heating  element,  which  is  located 
in  a 56  x 51  cm  plenum  chamber.  The  heating  capability  was 
not  necessary  for  this  work.  From  the  plenum  chamber  the 
flow  is  directed  via  a gradual  contraction  to  a smaller 
(15*2  x 30.4  cm)  plenum  chamber  with  flow  straightening 
elements.  The  flow  was  then  discharged  through  a 1.27  x 30.4 
cm  vertical  slot  on  a 91*^  x 30.4  cm  wall.  Twc  horizontal 
walls  (91-4  x 121.0  cm)  were  used  to  maintain  the  two  dimen- 
sionality of  the  flow. 

A two  dimensional  traversing  mechanism  was  used  to  move 
the  probe  apparatus.  The  range  of  the  traversing  device  was 
60  slot-widths  longitudinally  and  44  slot-widths  laterally, 
with  an  accuracy  better  than  ±0.04  cm.  The  measurements 
were  taken  at  x/D  = 35 > 4 5,  and  55  for  this  present  work. 

Mounted  on  the  probe  support  mount  was  an  L.  C.  Smith 
actuator  that  permitted  movement  of  the  sampling  probes  in 
the  lateral  direction.  The  detector  probes  were  movable 


only  by  moving  the  probe  support  mount.  The  actuator  pro- 
vided an  electrical,  remote  controlled  method  of  traversing 
the  sampling  probes.  A D.C.  voltage  across  a 10-turn  potent- 
iometer provided  a voltage  output  proportional  to  the  lateral 
traverse.  The  total  available  traverse  was  approximately 
30  slot-widths  relative  to  the  actuator  mount  with  an  accur- 
acy of  ±0.02  cm. 

The  probe  support  (see  Figures  11  and  12)  was  the  only 
part  of  the  experimental  setup  which  was  changed  from  that 
used  by  Jenkins  ( 197*0  • The  change  was  necessary  so  that 
two  sampling  probes  could  be  used.  The  two  probes  could  be 
moved  (in  the  lateral  direction,  manually)  up  to  a separation 
of  8 slot-widths.  The  sampling  probes  were  also  limited  to 
a maximum  separation  from  the  detector  probes  of  5 slot- 
widths.  Neither  of  these  constraints  were  found  to  be 
adversely  restrictive.  The  purpose  of  this  change  in  design 
was  to  provide  for  the  measurement  of  space-correlation 
coefficients  (see  Appendix  B). 

Measurements  were  taken  in  the  jet  to  determine  how 
well  it  agreed  with  work  from  previous  investigators.  The 
mean  velocity  and  relative  intensity  profiles  were  measured 
at  three  x/D  stations  and  are  shown  in  Figures  13  and  14 
respectively.  The  mean  velocity  profile  was  compared  with 
the  theoretical  results  of  Goertler  (from  Schlichting  ( 1968) ) 

1=1-  tanh2  (8) 

Um 

where  i\  = ay/x,  and  a is  a characteristic  of  the  jet.  For 
this  work  a was  found  to  be  9*28. 
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The  relative  turbulent  intensity  showed  agreement  with 
previous  work  by  other  investigators.  Bradbury  (1965)  showed 
a similar  plot  of  turbulent  intensity  increasing  from  the 
centerline  to  a maximum  value  at  y/b  = 0.7,  before  falling 
off  to  zero. 

The  turbulent  intensity  measured  at  the  centerline  of 
the  jet  was  compared  to  measurements  by  other  investigators 
(see  Figure  15)-  This  comparison  showed  these  measurements 
were  consistent  with  the  data  previously  reported. 

An  effort  was  made  to  assure  the  similarity  and  self- 
preservation  principles  of  the  jet  as  they  are  discussed  by 
Schlichting  (1968)  were  valid  in  the  case  of  this  setup,  one 
such  principle  states  that  the  half  width  of  the  velocity 
field  (y  = b at  T = ) is  directly  proportional  to  the  dis- 

tance downstream  of  the  jet  exit.  This  relationship  may  be 
stated  as 


b/D  = k x ( x/D  - Cx) 

where  k^is  the  widening  rate  of  the  velocity  field  and  C ^ is 
the  geometric  virtual  origin  of  the  jet. 

Figure  16  shows  the  half  width  widening  relationship  as 
measured  in  the  present  setup.  The  resulting  equation  is 

b/D  = 0.095U/D  + 0.789)  (9) 

Another  test  for  similarity  is  the  decay  of  the  center- 
line  velocity.  Theoretical  considerations  state  the  center- 
line  velocity  should  be  inversely  proportional  to  the  dis- 
tance downstream  of  the  jet  exit  squared.  This  relationship 


Variation  of  Turbulent  Intensity  Squared  Along  the  Centerline 
of  the  Jet 
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may  be  written  in  a linear  form,  as 

-2 

= k2(x/D  - C2) 

where  k2  is  the  decay  rate  of  the  centerline  velocity  and  C2 
is  the  kinematic  virtual  origin  of  the  jet.  As  noted  by 
Flora  and  Goldschmidt  (1969) > there  is  no  reason  to  expect 
any  relationship  between  the  two  virtual  origins  mentioned. 

Figure  17  shows  the  parabolic  decay  of  the  centerline 
velocity.  The  empirical  relationship  is  given  as 


0.l85(x/D  - 13.2) 


(10) 


The  empirical  constants  which  were  measured  in  this 
setup  are  compared  with  work  by  other  investigators  in  Table 

h. 

1.  For  Reynolds'  numbers  in  a range  of  1.0  to  8.0x10  , the 
variation  of  the  widening  rate  is  O.O85— 0.110  and  of  the  decay 
rate  of  the  centerline  velocity  is  0.150—0.36^. 

These  results  show  that  the  experimental  set  up  used  for 
this  work  is  in  general  agreement  with  previous  investigators. 
Other  tests,  especially  for  the  two  dimensionality  of  the 
flow  field  have  been  reported  previously  by  Jenkins  (197M 


and  were  not  repeated  by  this  author. 


2 . Procedure  to  Measure 

The  procedure  to  measure  the  lateral  velocity  of  the 
interface  is  described  below.  Figure  18  shows  a schematic 
of  the  data  measurement  procedure.  The  output  of  the  two 
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hot  wire  probes  which  are  located  in  known  positions  in  the 
flow  field  ((xQ,yo)  and  (xo,yQ-Ay)),  are  recorded  simultan- 
eously. Figure  2 shows  a typical  record  on  the  visicorder 
which  could  also  be  monitored  on  the  storage  oscilloscope. 
Comparing  the  two  simultaneous  traces,  the  length  of  time  it 
would  take  for  the  interface  to  travel  between  the  two 
probes  could  be  measured  as  the  time,  At,  between  the  onset 
of  turbulence  from  one  probe  to  the  onset  of  turbulence  at 
the  other  probe.  Averaging  over  sufficient  values  and  know- 
ing the  separation  of  the  two  probes,  Ay,  the  velocity  of 
the  interface  as  a function  of  the  probe  separation  and  the 
location  in  the  jet  may  be  determined  as 

v^x.y.Ay)  = 41 

The  positioning  of  the  probes  was  done  in  the  following 
manner.  Probe  2 (see  Figure  12)  was  positioned  at  the  x/D 
station  and  y/b  location  of  interest.  The  other  probe 
(probe  3)  was  then  positioned  at  different  lateral  separa- 
tions, Ay,  from  the  first  probe1.  Typically  a minimum  of 
three  values  of  Ay  were  used  to  find  V^(x,y,Ay).  From  them, 
the  extrapolated  V^, 

V; (x,y)  = Lim  V.(x,y,Ay) 

1 Ay-*-0  1 

was  computed.  This  corresponds  to  the  velocity  of  the  inter- 
face at  the  point  where  probe  2 was  located. 

1 Probes  2 and  3 mentioned  here  correspond  to  probes  2 and  1 
respecti\ely  of  Figure  1. 
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The  presence  of  folding  in  the  interface  is  a well  est- 
ablished fact.  Folding  is  encountered  when  a nonturbulent 
region  or  pocket  of  fluid  exists  closer  to  the  jet  axis  than 
another  turbulent  region  of  fluid.  One  example  of  an  inter- 
face configuration  which  would  lead  to  folding  is  shown  in 
Figure  19-  The  interface  configuration  is  meant  to  corre- 
spond to  the  hot  wire  traces  shown  in  the  figure. 

The  occurrence  of  folding  leads  to  difficulties  in  mod- 
eling the  interface,  as  the  interface  can  no  longer  be  as- 
sumed to  be  single  valued  (see  Paizis  and  Schwarz  (197*0)* 
The  percentage  of  folding,  ^ , is  a measure  of  the  occur- 
rence of  folding  in  the  interface.  It  is  the  fraction  of 
total  samples  observed  in  which  folding  is  noted,  and  is 
recognized  to  be  a function  of  location  in  the  flow  field. 

Another  problem  folding  introduces  is  an  ambiguity  in 
the  measure  of  the  velocity  of  the  interface.  From  Figure 
19  the  corresponding  different  values  of  would  be 

Vi1  = Ay/I^Ati  + + At5  + At6^  (H) 

Vi2  = a y/^At^  + At  g + At^  + At^  + At,,  + At^)  (12) 

- At 2 - At 2 + At^  + At^  + At^)  (13) 

is  a measure  of  the  velocity  of  the  interface  neglecting 
the  effects  of  folding.  V^2  obtains  a measure  of  the  velo- 
city of  the  interface  considering  the  effects  of  the  folding 
of  the  interface  in  magnitude  only.  considers  both  the 

direction  and  magnitude  of  the  folding  of  the  interface  in 


Figure  19  Possible  Form  of  Interface  With  Output  of  Hot 
Wire  Probes  Sh?wing  Folding 


its  measure  of  the  velocity  of  the  interface.  In  the  re- 
ported work,  all  three  cases  were  considered  although 
was  chosen  to  be  the  most  meaningful  measure  of  the  velocity 
of  the  interface. 
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III  EXPERIMENTAL  RESULTS 


1 . The  Velocity  of  the  Interface 

The  measurement  of  the  velocity  of  the  interface  has 
been  carried  out  as  per  the  procedure  given  in  the  previous 
section.  The  results  are  shown  and  explained  in  this  sect- 
ion. 

Figure  20  shows  the  mean  velocity  of  the  interface,  V^, 
(specifically  °f  Equation  (l3))asa  function  of  the  lat- 
eral displacement  in  the  plane  jet  for  three  different  x/D 
stations.  This  figure  shows  a monotonic  decrease  of  as 
it  moves  away  from  the  centerline  of  the  jet  (y/b=0).  It 
should  also  be  noted  that  no  strong  variation  is  seen  with 
changes  in  axial  location.  This  could  be  due  to  the  limited 
range  of  x/D  values  or  due  to  experimental  scatter.  In  real- 
ity, similarity  suggests  that  should  scale  with  and/or 

bfYm’ 

There  is  some  question  as  to  the  best  method  of  nondi- 
mens ionalizing  V ^ . Figures  21,  22  and  23  show  plotted  in  non- 
dimensional  terms.  In  the  first  plot,  is  made  nondimen- 
sional  by  b,  the  half  width  of  the  velocity  field  and  fym# 
the  maximum  crossing  frequency.  This  grouping  is  suggested 
by  the  analysis  used  to  predict  the  value  of  the  velocity  of 

the  interface.  The  second  plot  shows  V.  scaled  by  TL,  the 

x m 


mean  centerline  velocity.  This  grouping  seems  to  give  a 
reasonable  scaling  of  the  data.  It  shows  the  order  of  mag- 
nitude of  to  range  from  0.8  to  0.1  of  the  mean  centerline 
velocity  while  moving  through  the  intermittent  region  of  the 
jet.  The  third  plot  shows  made  nondimens ional  by  0,  the 
local  mean  velocity.  From  this  figure  we  see  that  the  ratio 
cf  Vj/ff  increases  as  one  moves  away  from  the  centerline  of 
the  jet. 

Figures  24,  25  and  26  show  the  same  information  as  a 
function  of  intermittency . From  Figure  25  it  can  be  seen 
that  a great  deal  of  the  variation  in  the  ratio,  V^/TJm, 
occurs  near  the  edges  of  the  intermittent  region  of  the  jet. 

Figures  2?  and  28  compare  the  experimental  values  of 
with  that  predicted  from  the  analysis.  The  former  shows  the 
information  as  a function  of  y/b  for  different  x/D  stations, 
while  the  latter  is  shown  versus  intermittency. 

It  seems  that  the  measured  values  of  V\  (chosen  as  Equa- 
tion (13))  do  agree  in  both  magnitude  and  in  general  trend 
with  the  analysis.  The  variation  occurs  at  the  edges  of  the 
intermittent  region  where  measurements  are  more  difficult  to 
take. 

The  two  other  forms  of  V\ , which  were  defined  by  Equa- 
tions (11 ) and  (12)  and  required  due  to  the  existence  of 
foldover,  are  compared  with  the  values  predicted  from  the 
analysis  in  Figures  29  and  30 . Figure  29  shows  in  non- 

dimensional  terms  compared  with  the  predicted  value  of  V^. 
may  be  thought  of  as  the  velocity  of  the  interface 
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without  considering  the  folding  of  the  interface.  As  shown 
in  Figure  29,  the  measured  value  of  appears  to  be  half 
the  predicted  value  of  V\. 

In  a similar  manner  Figure  30  compares  V^2  with  the 
predicted  value  of  V.  V.„  is  the  velocity  formed  by  taking 
the  absolute  value  of  each  time  increment,  At.  Thus  the 
contribution  due  to  foldover  is  positive  rather  than  negative 
as  it  is  for  (or  from  Equation  ( 13) ) • The  difference 
between  the  values  of  and  V^2  is  a measure  of  the  asym- 
metric nature  of  the  interface  and  its  folding. 

2 . Further  Measurements  of  the  Interface 
In  addition  to  the  measurement  of  the  velocity  of  the 
interface,  other  results  were  obtained  which  will  be  pre- 
sented here. 

a.  Fronts  and  Backs  of  the  Interface 
In  the  literature  concerning  the  turbulent-nonturbulent 
interface,  the  terms  front  and  back  of  the  interface  are 
used.  The  front  of  the  interface  is  defined  as  the  part  of 
the  interface  where  the  flow  goes  from  nonturbulent  to  turb- 
ulent. This  view  comes  from  the  idea  of  the  interface  which 
separates  turbulent  and  nonturbulent  fluid  being  convected 
downstream  by  the  mean  flow.  Similarly  the  back  of  the  in- 
terface is  the  part  where  the  flow  goes  from  turbulent  to 
nonturbulent  fluid.  This  may  be  seen  in  Figure  19  which 
shows  a possible  shape  of  the  interface. 
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If  the  model  of  the  interface  is  taken  as  that  of  the 
interface  moving  laterally  in  a varying  periodic  motion, 
then  the  front  of  the  interface  would  correspond  to  the  part 
of  the  interface  indicated  by  the  time  increments  and 

at*,  in  Figure  19*  This  is  an  apparent  outward  motion  (away 
from  the  jet  centerline)  of  the  interface.  In  a similar  man- 
ner, the  back  of  the  interface  would  correspond  to  the  time 
increments  at^  and  atg  of  Figure  19  and  could  be  interpreted 
as  an  inward  motion  of  the  interface. 

The  point  of  discussion  is  that  the  fronts  and  backs  of 
the  interface  are  not  necessarily  identical  (symmetric). 
Paizis  and  Schwarz  (197*0  discuss  this  for  their  work  in  a 
wall  jet  and  have  found  that  the  slope  of  the  front  of  the 
interface  is  steeper  than  the  slope  on  the  back  of  the  inter- 
face. Oswald  and  Kibens  (1971)  have  also  found  the  fronts 
and  backs  of  the  turbulent-nonturbulent  interface  to  differ 
in  the  wake  of  a circular  disk.  They  have  taken  measurements 
which  show  that  the  convective  velocity  of  the  front  of  the 
interface  is  faster  than  the  convective  velocity  of  the  back. 
Figure  31  compares  the  velocity  of  the  front  of  the  inter- 
face with  the  velocity  of  the  back.  This  shows  that  the 
interface  appears  to  move  outward  faster  than  it  moves 
inward.  This  could  be  a result  of  the  difference  in  the 
slopes  of  the  front  and  back  of  the  interface  which  is  dis- 
cussed in  the  work  of  Paizis  and  Schwarz  (1974). 
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b.  Folding 

The  possible  presence  of  folding  in  the  interface  has 
been  a source  of  uncertainty  in  conditional  measurements. 

Its  presence  in  the  case  of  a wall  jet  was  first  noted  by 
Paizis  and  Schwarz  (197^).  although  ambiguity  in  the  gate 
level  of  their  turbulence  detectors  could  have  lead  to  the 
apparent  presence  of  folding.  In  this  work  the  observation 
of  simultaneous  traces  such  as  shown  in  Figure  2 has  convinc- 
ingly indicated  the  presence  of  folding  (as  per  the  schemat- 
ic of  Figure  19) • 

The  parameter,  was  defined  as  the  percentage  of 

folding  in  the  interface.  ^ is  shown  in  Figure  32  as  a 
function  of  the  location  in  the  jet.  is  also  shown  in 

Figure  33  for  both  the  front  and  back  of  the  interface  as  a 
function  of  location  in  the  jet.  From  these  two  figures  it 
seems  that  folding  is  more  likely  to  occur  toward  the  center- 
line  of  the  jet  (in  the  intermittent  region)  and  more  like- 
ly to  be  detected  on  the  front  of  the  interface.  It  should 
be  noted  that  there  can  be  no  folding  outside  the  intermit- 
tent region. 
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Figure  33  Percentage  Folding  for  Fronts  and  Backs  of 
the  Interface  vs.  y/b 
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IV  DISCUSSION  OF  RESULTS 


This  work  has  been  mostly  experimental  with  the  empha- 
sis placed  on  measurements  taken  in  the  laboratory.  However 
there  is  one  analytical  result  which  has  come  from  this  work 
and  deserves  to  be  mentioned  here.  The  expression  for  the 
velocity  of  the  interface, 

-2f„ 


Vi  W/l y 


(6) 


which  was  derived  in  Chapter  II  from  basic  parameters  of 
turbulent  flows  was  compared  for  different  shear  flows. 

There  was  a trend  of  similarity  noted  for  the  flows  analyzed. 
This  is  significant  in  that  this  is  the  first  time  a charac- 
teristic velocity  of  the  interface  has  been  identified  show- 
ing similarity  for  different  situations. 


1 . Experimental  Accuracy 

In  order  to  quantify  the  accuracy  of  the  measurements 
reported,  an  analysis  of  the  experimental  accuracy  of  the 
measurements  will  be  presented. 

Applying  an  uncertainty  analysis  to  or  V^(Ay),  the 
primary  measured  quantity,  an  estimate  of  the  uncertainty 
may  be  obtained  (Appendix  C shows  the  method  used  to  obtain 
from  the  data  taken  from  the  visicorder  traces)*  the 


source  of  error  in  the  measurements  to  obtain  V^(Ay)  comes 
from  uncertainty  in  the  measure  of  At  from  the  visicorder 
charts  and  Ay,  the  probe  separation.  For  example,  the  un- 
certainty in  V^(Ay= 0.254)  based  on  a single  measure  of  At 
would  be  ± 37.2#  (this  is  the  sum  of  6At/At  = ±29.3$  and 
6 Ay/ Ay  = ±7.9$).  If  all  the  time  increments,  At,  were  of 
the  same  duration,  the  uncertainty  due  to  the  measure  of  At 
would  be  reduced  by  a factor  of  the  inverse  of  the  square 
root  of  the  number  of  samples  taken.  Using  this  factor 
the  uncertainty  of  V^(Ay=0.254)  is  reduced  to  ±12.2#  and 
the  uncertainty  of  V^(x,y)  found  by  extrapolation  is  ±13. 5#. 

The  fact  is  that  the  time  increments  are  not  all  of  the 
same  duration.  Hence,  in  order  to  obtain  a meaningful  aver- 
age, a large  number  of  samples  must  be  recorded.  The  actual 
number  of  samples  recorded  was  limited  by  two  constraints. 

One  the  laborious  manual  process  employed  in  collecting  the 
samples,  and  secondly  their  distribution  due  to  the  physical 
location  of  the  probes  in  the  intermittent  region. 

If  it  is  assumed  that  the  measured  values  of  At  follow 
a normal  distribution,  then  an  estimate  of  the  further  error 
in  V^(x,y)  due  to  the  number  of  samples  taken  to  obtain  the 
mean  value  can  be  obtained. 

Data  was  taker,  at  eleven  points  throughout  the  fully 
developed  region  of  the  jet  and  at  a total  of  36  settings  of 
(x,y,Ay).  The  mean  number  of  samples  of  the  time  increment, 
At,  taken  at  each  setting  was  38  with  a range  of  12  to  74 
samples.  The  number  of  samples  taken  was  dependent  primarily 
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on  two  factors:  the  y/b  location  in  the  jet  and  the  probe 
separation.  The  determination  of  At  from  the  visicorder 
traces  was  difficult  for  large  probe  separations  and  at  the 
edges  of  the  intermittent  region  (y^O.O  or  y=1.0).  It 
should  be  noted  that  at  the  smallest  probe  separations  the 
number  of  samples  was  consistently  large. 

In  an  effort  to  determine  how  the  number  of  samples 
affected  the  accuracy  of  At  ( and  hence,  Vt(x,y)),  the  stan- 
dard deviation  of  the  mean,  am , given  as 

^ = VT 

was  calculated,  a is  the  standard  deviation  of  the  set  of 
samples  measured  and  n is  the  number  of  samples  in  the  set. 
For  a typical  case  the  ratio  am/V t was  found  to  be  O.J. 

This  then  implies  a standard  deviation  of  the  mean  value  of 
V^(ay)  due  to  the  number  of  samples  taken  to  be  ±30#  and  a 
corresponding  probable  error  of  about  20#.  This  error  is  in 
addition  to  the  uncertainty  calculated  above. 

From  this  analysis  the  measurements  of  Vt(x,y)  taken 
and  reported  in  this  work  show  an  uncertainty  of  ±13-5^  due 
to  the  measurement  procedure  and  a probable  error  of  ±20# 
due  to  the  number  of  samples  taken.  This  would  give  an  er- 
ror band  of  ±33#  about  the  measured  values  of  V^(x,y).  This 
seems  quite  reasonable  for  the  type  of  measurement  taken. 


2 . The  Asymmetric  Nature  of  the  Interface 
One  characteristic  of  the  interface  which  has  been 


brought  forth  by  this  work  is  its  lack  of  symmetry.  This  was 
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mentioned  previously  when  V^,  V^2,  and  were  introduced 
(see  Equations  (11),  (12),  and  (13)  and  Figures  28  through 
30).  The  difference  between  and  is  expected  since 

V.^  accounts  for  the  negative  components  of  At  which  ne- 
glects. However  the  difference  between  V^2  and  might 
no1  be  exj'Pcted  since  V^2  is  simply  the  magnitude  of  the 
velocity  without  regard  to  direction.  If  the  magnitude  of 
At  due  to  folding  was  the  same  as  the  normal  values  of  At, 
then  would  be  the  same  as  V.^.  This  is  not  the  case  as 
can  be  seen  from  a comparison  of  Figures  29  and  30.  This 
shows  that  the  velocity  of  the  interface  during  folding  is 

V 

different  than  that  without  folding. 

Another  aspect  of  the  asymmetric  nature  of  the  inter- 
face is  the  marked  difference  between  the  velocities  of  the 

I 

front  and  back  of  the  interface.  Paizis  and  Schwarz  (1974) 
in  their  work  in  a turbulent  wall  jet,  noted  that  the  slope 
of  the  front  of  the  interface  (dy^/dx)  was  greater  than  the 
slope  of  the  back  of  the  interface.  This  is  in  agreement 
with  this  current  work.  The  velocities  of  the  interface  for 
the  front  and  back  of  the  interface  (see  Figure  31)  may  be 
related  to  the  slope  of  the  interface  by  the  relationship 

ayj_  dy-j/dt  v. 

dx  ” dx/9 1 U 

This  has  been  done  for  x/D  = 35  and  is  shown  in  Figure  3^* 

The  figure  clearly  shows  the  difference  between  the  front 
and  back  of  the  interface  when  considering  an  average  shape 
including  all  folding  and  non-folding  cases. 

i 

I 


Average  Shape  of  the  Interface 


60 


Most  of  the  difference  occurs  at  Y^O.5  which  coincides  with 
the  region  of  the  jet  that  exhibits  the  most  folding  in  the 
interface . 


3 • Folding  and  Entrainment 

One  of  the  secondary  goals  of  this  work  was  to  provide 
more  information  about  the  entrainment  of  ambient  fluid  into 
the  turbulent  flow. 

It  has  been  proposed  by  Roshko  and  discussed  by  Davies 
and  Yule  (197*0  that  the  entrainment  process  involves  engulf- 
ment  of  ambient  fluid,  particularly  by  coalescing  eddies. 

This  idea  leads  to  the  conclusion  that  entrainment  would  be 
dependent  upon  the  motion  of  the  interface  rather  than  vis- 
cous mechanisms. 

A quantification  of  the  motion  of  the  interface  and 
especially  the  effect  of  folding  in  the  interface  would  seem 
to  be  important.  One  question  which  has  been  raised  by  this 
work  and  whose  answer  might  give  more  information  about  en- 
trainment is  why  does  the  detection  of  folding  increase  as 
one  moves  toward  the  centerline  of  the  jet?  At  present  this 
can  only  be  stated  as  an  experimental  observation. 

Conclusions 

Overall  this  work  seems  to  have  provided  a unique  mea- 
sure of  the  velocity  of  the  interface  in  the  lateral  direc- 
tion. This  was  accomplished  by  viewing  the  motion  of  the 
interface  in  a manner  different  from  earlier  investigators. 
Other  characteristics  of  the  interface  have  also  been 


62 


V FUTURE  WORK 

It  was  felt  that  there  should  be  some  discussion  of  the 
future  work  to  be  undertaken  as  a continuation  of  this  cur- 
rent work.  There  are  two  areas  which  it  is  felt  are  of  in- 
terest and  will  be  investigated  in  the  future. 

The  first  area  is  the  interface  and  its  characteristics. 
There  has  been  no  real  attempt  made  here  to  model  the  inter- 
face, although  the  view  of  the  interface  given  by  may 
lead  naturally  to  a model.  There  is  also  a need  to  apply 
some  model  to  the  problem  of  entrainment.  This  should  be  one 
of  the  goals  in  any  future  work  in  addition  to  more  experi- 
mental measurements  of  the  interface. 

The  second  area  of  interest  is  concerned  with  condition- 
al averaged  measurements  in  a two  dimensional  plane  jet.  As 
an  extension  of  the  work  by  Jenkins  (197^),  the  intention  is 
to  take  measurements  of  the  lateral  scales  of  turbulence  in 
the  fully  turbulent  region.  The  result  is  a measure  of  the 
scales  with  respect  to  the  interface  at  different  locations 
in  the  jet. 

The  purpose  of  these  latter  measurements  is  to  shed 
some  insight  on  free  shear  flows  and  answer  the  question — 
how  does  shear  change  the  turbulent  structure  within  the 
turbulent  region  itself?  Appendix  B contains  more  informa- 
tion on  these  measurements. 
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APPENDIX  A THE  DESIGN  OF  A LOW  SPEED  WIND  TUNNEL 


This  appendix  is  meant  to  document  the  work  done  by  the 
author  in  the  design  and  construction  of  a low  speed  wind 
tunnel.  The  purpose  of  the  wind  tunnel  is  to  be  part  of  an 
experimental  setup  for  taking  mean  and  conditional  measure- 
ments of  temperature  and  velocity  in  the  cold  wake  of  a cyl- 
inder . 

The  final  design  of  the  wind  tunnel's  test  section  is 
shown  in  Figure  A-l , and  Figures  A-2,  A~3  and  A-4  show  the 
actual  wind  tunnel  built  by  the  shop  staff  of  the  Herrick 
Laboratories . 

The  test  section  has  a cross  section  of  30-48  x 45-72  cm 
wide  and  2.44  m in  length.  One  side  of  the  test  section  is 
made  of  1.91  cm  plexiglas  with  a 1.27  cm  slot  at  the  mid- 
plane of  the  tunnel.  The  slot  is  to  be  used  to  insert  probes 
into  the  flow  field. 

The  tunnel  was  designed  as  an  open  loop  wind  tunnel 
with  the  blower  at  the  exit,  sucking  the  air  through  the  tun- 
nel. This  was  done  to  achieve  the  best  results  with  a mini- 
mum of  cost  and  labor. 

A tunnel  such  as  the  one  built  will  have  a relatively 
low  turbulent  intensity  level  in  the  test  section  without 
the  need  for  screens,  straws  or  honeycombs.  This  was  the 
major  reason  for  the  choice. 


Figure  A-l  Schematic  of  Wind  Tunnel 
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The  tunnel  is  equipped  with  a LAU  blower,  model  A10-10 
AC,  squirrel  cage  type.  It  is  powered  by  a 1.12  kW  motor 
rated  at  3450  RPM. 

In  the  design  stages,  a flow  rate  of  1.7  m^/sec  (3600 
cfm)  was  chosen  as  a specification  so  that  the  tunnel  would 
have  a free  stream  velocity  of  9.14  m/sec  (30  fps).  This 
was  the  minimum  velocity  which  would  be  usable.  After  the 
construction  of  the  tunnel  was  completed,  the  maximum  free 
stream  velocity  was  measured  with  a pitot  probe  and  was 
found  to  be  10  m/sec.  The  free  stream  velocity  can  be  con- 
trolled by  a flapper  valve  that  covers  the  exhaust  of  the 
blower. 

Since  the  wind  tunnel  was  meant  to  be  used  to  take 
measurements  in  the  wake  of  a cooled  circular  cylinder,  a 
cooling  system  was  designed  and  partially  built  for  the  tun- 
nel. Figure  A-5  shows  a design  for  the  cooling  system.  It 
is  proposed  that  Freon  502  (BP  = -45°  C ) be  circulated 
through  the  system  and  hence  maintain  the  cylinder  within 
the  tunnel  at  a constant  temperature  of  -45°  C. 

Calculations 

Calculations  were  made  prior  to  the  construction  of  the 
wind  tunnel  to  make  sure  the  desired  measurements  would  be 
theoretically  possible  in  the  setup  as  designed.  The  calcu- 
lations will  be  reviewed  here. 

To  establish  Reynolds  number  similarity,  Townsend  (1956) 
has  found  practically  no  effect  on  the  mean  velocity  profiles 
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when  Red  > 800.  The  Reynolds  number  based  on  the  cylinder 
diameter  was  calculated  with  tJo  = 9.1^  m/sec  and  d = 0.635  cm 
and  was  found  to  be  Re^  = 3800  for  an  air  temperature  of  23° 

C (v  = 15.2  x 10  m /sec).  It  should  be  noted  that  if  the 
effective  temperature  is  lowered,  then  the  Reynolds  number 
will  be  increased. 

Another  question  to  answer  was  an  estimate  of  the  widen- 
ing rates  of  both  the  velocity  and  temperature  profiles,  and 
also  an  estimate  of  the  velocity  and  temperature  defect  mag- 
nitudes at  different  stations  downstream. 

It  is  known  that  the  half  width  of  the  velocity  field, 

b , (where  db^/dx  is  the  widening  rate)  is  proportional  to 

a 

the  square  root  of  the  distance  downstream,  x2 , in  the  simi- 
larity region  (see  Schlichting  (1968) ) . Thus  the  form  of  the 
half  width  may  be  written 

/x  - x \ i 

Vd  = ci( — n f*-1’ 

where  xq  is  the  virtual  origin  of  the  wake. 

Using  Townsend's  (19^9)  data  discussed  in  Hinze  (1959). 
the  velocity  defect  profile  in  a wake  may  be  approximated  by 

°l/0lm  = exP(-(orfe,2)  (A-2> 

where  a nond imensional  grouping  of  y.  From  this  equa- 

tion the  half  width  of  the  velocity  field  is  found  to  be 

bu/d  = o.213(2L-^-§)^  ( A— 3 ) 

Making  the  assumption  that  the  same  virtual  origin  will 
exist  in  the  designed  test  setup  as  was  found  by  Townsend  in 
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his  experimental  work  (a/d  = 20),  it  is  possible  to  estimate 
the  half  width  of  the  wake  as  a function  of  the  distance 
downstream.  The  widening  rate  of  the  velocity  field  is 
shown  in  Figure  A-6. 

In  a similar  fashion  the  half  width  of  the  temperature 
field  in  the  wake  cf  a circular  cylinder  will  have  the  form 


b0/d  = C 


(A-4) 


From  the  temperature  measurements  taken  by  Townsend 
(1949)  and  analyzed  by  Hinze  (1959).  it  has  been  found  that 
a reasonable  fit  to  the  data  is  obtained  from  a model  based 
on  a constant  coefficient  of  thermal  transport.  From  this 
model  the  temperature  profile  is  a Gaussian  solution  similar 
to  the  velocity  profile, 

o/%  = exP<-<oTFS>2>  (A’5> 

The  widening  rate  of  the  temperature  field  may  then  be 
found  to  be 


b0/d  = 0.290 


i 


(A-6) 


By  making  the  same  assumptions  concerning  the  virtual 
origins  of  the  wake,  an  estimate  of  the  spreading  rate  of 
the  temperature  may  be  found.  The  widening  rate  of  the  temp- 
erature profile  is  shown  in  Figure  A-6.  The  two  widening 
rates  are  shown  together  for  the  sake  of  comparison.  The 
figure  shows  the  temperature  spreading  faster  that  the  mean 
velocity.  This  difference  in  spreading  rates  has  been 
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explained  by  Jenkins  (197^)  for  a plane  jet  to  be  caused  by 
the  different  slopes  of  the  point  averaged  velocity  and  temp- 

^ A 

erature  profiles  (de/dy  and  dfl/dy)  at  the  interface. 

It  is  also  desired  to  estimate  the  velocity  and  tempera- 
ture defects  on  the  centerline  of  the  wake  as  a function  of 
the  distance  downstream.  According  to  Schlichting  (1968),  the 
centerline  velocity  defect  varies  with  the  distance  down- 
stream as 

Uln/Uo  = ’‘Hr*)’4  <A-?> 

in  the  similarity  region.  From  data  taken  by  Townsend  (1956) 
the  velocity  defect  relationship  may  be  approximated  as 

Uo/Ulm  = 1 • 09(x/d  - 25)'  ( A-8 ) 

Since  no  data  was  readily  available,  the  temperature 
defect  at  the  centerline,  0^,  was  approximated  by  forming  an 
energy  balance.  The  heat  out  of  the  system  is  found  by  the 
heat  transfer  to  an  infinite  cylinder  in  a crossflow  situa- 
tion. The  Reynolds  number  based  on  the  diameter  of  the  cyl- 
inder and  the  bulk  temperature  of  the  fluid  is  5300.  From 
Holman  (1968)  the  Nusselt  number  relation  corresponding  to 
the  Reynolds  number  range  is  given  as 

Nud  = 0.17MRed)0,618  (A-9) 

For  the  given  flow  situation  the  Nusselt  number  is  38.1, 
which  corresponds  to  a heat  flux  per  unit  length  of 

Qout  = 167*3  (W/m)  (A-10) 
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The  heat  flux  into  the  system  comes  from  the  entrain- 
ment of  ambient  fluid  into  the  wake.  Across  a slice  of  the 
flow,  the  heat  flux  is  given  by 

Qin  = / fUcr)0dy  (A-ll) 

where  U is  the  velocity  of  the  fluid  in  the  x direction  and 
0 is  the  temperature  defect. 

By  using  the  definition  of  velocity  defect,  = UQ  - U, 
and  assuming  the  velocity  and  temperature  defect  profiles 
may  be  described  as  Gaussian  solutions, 

U1  = ulmexP("a2y2)  (A-12a) 

0 = 0mexp(-32y2)  ( A-12b) 

with  a = -Jin  2/bu  and  g = -Jin  2/bQ , we  may  solve  for  e . 

Combining  these  equations,  the  heat  flux  into  the  sys- 
tem is 


f°pVmLe’B  y dy  - fcpUlr>9m'^e'(°l2+e2)y2dy  U'13) 

The  solution  of  (A-13)  gives 


Q in  = Pc~u~9. 


Qin  = fcpV"(^  ’ "r)  <A-l4) 

2 2 A 
where  y = (a  + 0 )s. 

The  heat  flux  in  and  out  of  the  system  may  be  equated, 
since  the  flow  is  assumed  steady  and  there  are  no  losses. 
Equating  (A-10)  and  (A-14),  gives  an  expression  containing 

62^1 
U 

o 

0 
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/U  ui  \-l 

% = °-°783{f  - -iS)  (A-15) 

The  temperature  defect  along  the  centerline  of  the 
wake  may  be  estimated  from  (A-15)*  This  is  shown  in  Figure 
A- 7 along  with  the  velocity  defect  along  the  centerline  of 
the  wake . 

Conclusions 

A sample  of  some  of  the  results  that  may  be  obtained 
from  the  above  calculations  are  found  in  Table  A-l.  Case  I 
of  the  table  was  taken  as  a standard  set  of  variables  which 
might  be  used  in  the  wind  tunnel.  Case  II  shows  what  changes 
might  be  expected  as  a result  of  changing  the  cylinder  dia- 
meter from  0.635  cm  to  O.3I8  cm.  Case  III  shows  what  changes 
might  occur  with  a change  of  cooling  agents  (from  Freon  502, 
BP  - -45°  C to  liquid  N,)(  BP  = -195.8°  C). 

The  information  which  can  be  gained  from  these  calcula- 
tions, given  a set  of  variables  is: 

1.  What  is  the  maximum  x/d  station  at  which  measure- 
ments may  be  taken? 

2.  At  x/d„,ov,  what  is  the  width  of  the  velocity  and 

m 3.x 

temperature  field?  Will  there  be  any  interference  from 
the  sides  of  the  wind  tunnel? 

3.  At  x/d__  . what  is  the  magnitude  of  the  velocity 
defect  and  the  temperature  difference?  Can  these  mag- 
nitudes be  measured  accurately? 
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Table  A-l  Predicted  Results  for  Wind  Tunnel 
Control 


Variables 

Case  I 

Case  II 

Case  III 

Uq  (m/sec) 

9.14 

9.14 

9.14 

CD 

o 

o 

o 

67.8 

67.8 

218.3 

d (cm) 

0.635 

0.318 

0.635 

Parameters 

Red 

5300 

2700 

8800 

x/Vx 

360 

700 

360 

xmax  <m) 

2.29 

2.23 

2.29 

Predicted 

Measurements 

@ x/d 
max 


b^  (cm) 

2.64 

1.82 

2.64 

bfl  (cm) 

3-59 

2.47 

3.59 

Ulm  (m/sec) 

0.458 

0.323 

0.458 

<°  C> 

0.205 

0.294 

0.660 

L t 
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APPENDIX  B CONDITIONAL  MEASUREMENT  OF  A 

This  appendix  is  intended  to  present  the  preliminary 
results  of  an  attempt  to  measure  the  point  averaged  lateral 
space  scales  in  a two  dimensional  plane  jet.  It  is  intended 
to  discuss  the  circuitry  necessary  to  measure  the  point  aver- 
aged space  correlation,  to  show  the  preliminary  results  ob- 
tained at  one  location  in  the  jet,  and  to  discuss  the  signi- 
ficance of  the  proposed  measurements. 

A point  averaged  measurement  is  a special  type  of  cond- 
itional measurement.  When  the  interface  is  sensed  to  be  in 
a predetermined  position  a sample  from  another  probe  is  taken. 
The  average  of  a large  number  of  these  samples  is  a point 
averaged  value.  It  is  recognized  as  the  mean  value  in  the 
fully  turbulent  region,  at  a specified  distance  from  the 
interface,  when  the  interface  is  in  a predetermined  position. 

Measurement  Procedure 

Figure  B-l  shows  a schematic  of  the  procedure  used  for 
measuring  the  mean  space  correlation,  g(r),  and  the  point 
averaged  correlation  coefficient.  It  shows  the  signals  taken 
from  two  hot  wire  probes  (see  Figures  11  and  12  for  position- 
ing of  the  probes),  which  are  separated  by  a distance,  r,  in 
the  lateral  direction.  These  signals  are  amplified,  AC 


Figure  B-l  Schematic  of  A Measurement  Process 
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coupled,  and  finally  multiplied  together  using  an  Analog 
Devices  multiplier  circuit,  model  426  K.  The  conditioned 
signal  is  then  biased  by  a known  constant  DC  voltage  and 
measured  using  a B & K RMS  meter  and  CALICO  digital  voltmeter. 

The  lateral  space  correltation  coefficient,  g(r),  is  de- 
fined by  Taylor  (1935)  as 


u(?1)u(S1+r) 

g(r)  = 


( B— 1 ) 


This  is  equivalent  to 


g<r)  = i“?£7 

e3  e4 


(B-2  ) 


where  e^  and  e^  are  the  voltage  fluctuations  corresponding 
to  the  velocity  fluctuations,  u(?1)  and  u(T^+r)  respectively. 
Thus  from  the  schematic  the  correlation  coefficient  may  be 

3.  Id 

obtained  from  the  measurement  of  ae^' , be^’  and  fo®3^4  as 

ab- 

(B-3) 


•<*>  ■ 4-^t) 


The  lateral  space  macroscale, A , is  defined  as 

o 


Ag  =Jo  S(r5dr 


(B-4) 


As  a means  to  establish  a precedent  the  macroscale  in  this 
work  will  be  measured  as 


Ag  =/or'e<r)dr 

where  g(r’)  = 0,  since  the  possiblity  exists  for  g(r)  to  be 
a negative  value.  This  relieves  the  ambiguity  of  the  defin- 
ition (B-4). 
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The  signal  from  the  voltage  summer,  Vdc  + is  AC 

coupled  to  eliminate  the  bias  voltage  and  put  into  the  vari- 
able time  delay  before  going  to  the  Waveform  Eductor.  The 
purpose  of  ihtroducing  a delay  time  is  to  match  the  delay 
time  introduced  by  the  interface  detector  circuit. 

The  signal  I(t)  which  triggers  the  Waveform  Eductor 
comes  from  the  interface  detector  circuit  and  the  two  hot 
wire  detector  probes.  The  output  of  the  detector  is  a sig- 
nal, I(t),  which  is  unity  when  the  flow  seen  by  the  detector 
probes  is  turbulent  and  zero  when  the  flow  is  nonturbulent . 
The  time  derivative  of  I(t)  will  give  a nonzero  reading  at 
the  times  when  the  interface  is  crossing  the  detector  probes. 
The  intricacies  of  the  interface  detector  circuit  are  dis- 
cussed thoroughly  by  Jenkins  (1974). 

When  the  Waveform  Eductor  triggers  on  the  signal  I(t), 
a sample  of  the  signal  (in  this  case  ^e^e^)  'taken*  Af- 
ter an  appropriate  number  of  samples  is  taken  and  averaged, 

3.  Id  3.  Id 

the  result  is  Yoe3e4»  the  Point  averaged  value  of  Yoe3e4‘ 

The  point  averaged  correlation  coefficient,  £(r),  is 
defined  as 

t(r)  = (B-6) 

e3  e4 

where  e^'  and  e^'  are  the  root  mean  squared  values  of  the 

A 

values  e^  and  (ie.  from  £ *»  E + e).  The  circuitry  to 
make  these  measurements  has  not  been  assembled,  although  it 
should  not  be  difficult  to  do  using  digital  techniques. 
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In  order  to  approximate  the  point  averaged  lateral  space 
scale,  an  approximate  point  averaged  correlation  coefficient, 


g1(r),  is  defined  as 


= 


e364 


e3e^max 


(B-7) 


and  then  the  approximate  point  averaged  lateral  scale  will 


be  given  by 


Ki = £'Mr)dr 


(B-8) 


To  form  an  idea  of  the  validity  of  this  approximation, 
the  same  approximation  is  made  with  regard  to  the  mean  cor- 
relation coefficient*  That  is 


K (r)  = -55- 

1 ^max 


(B-9) 


will  be  used  to  define  an  approximate  lateral  space  macro- 


scale , 


-V  = /0r'®i(r,dr 


(B-10) 


which  will  be  compared  to  Ag*  This  will  give  some  basis  to 
judge  the  accuracy  of  the  approximation,  Agl* 

Preliminary  Results 

There  is  some  preliminary  data  to  report.  The  measure- 
ments described  above  were  made  at  position  in  a two  dimen- 
sional plane  jet  (x/D  = 35  and  y/b  = 1.65)* 

Figure  B-2  shows  a comparison  of  the  mean  lateral  macro- 
scale and  the  approximate  macroscale.  Although  these  are 
preliminary  results,  it  does  not  seem  thatA^  agrees  very 
well  with  A • 

o 
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APPENDIX  C EXAMPLE  OF  PROCEDURE  USED 
TO  DETERMINE  Vi  FROM  DATA 

This  appendix  will  show  the  steps  taken  to  determine 
from  the  raw  data  taken  from  the  visicorder  traces.  As  an 
example  the  actual  data  for  one  point  in  the  flow  field, 
x/D  = 35  and  y/h  = 2.25,  will  be  analyzed.  Table  C-2  shows 
the  raw  data  as  it  was  taken  from  the  visicorder  traces. 

The  average  time  increments,  a^  and  are  found  as 

the  arithmetic  mean  of  the  samples  and  At  is  given  as 

At  = J(Atf  + Ztb) 

The  corresponding  velocities,  V^,  V^.  and  Vi  can  -fc^en  136 
computed  as  a function  of  the  probe  separation.  These  val- 
ues are  shown  in  Table  C-l. 

Table  C-l  Measured  Values  of  Zt  and  V(Ay)  at  x/D  = 35  and 
y/b  = 2.25 


ay  (mm) 

2.5  . _ 

5.1 

7.6 

(msec) 

1.0? 

3.20 

3.60 

2.25 

3-53 

3.78 

1.66 

3-36 

3.69 

Vfi  (m/s) 

2.37 

1.58 

2.12 

vbi 

1.13 

1.44 

2.01 

1.53 

1.51 

2.07 
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